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Phonon-assisted singlet-singlet relaxation in semicon- 
ductor quantum dot molecules is studied theoretically. 
Laterally coupled quantum dot structures doped with two 
electrons are considered. We take into account interac- 
tion with acoustic phonon modes via deformation poten- 
tial and piezoelectric coupling. We show that piezoelec- 
tric mechanism for the considered system is of great im- 
portance and for some ranges of quantum dot molecule 
parameters is the dominant contribution to relaxation. 



1 Introduction Quantum dot molecules (QDMs) have 
attracted much interest in experimental and theoret- 
ical ID investigations due to their proposed application in 
various implementations of quantum information process- 
ing schemes. In particular, two-electron spin states in cou- 
pled quantum dots (QDs) have been employed in a num- 
ber of realizations of quantum gates 1|4,5 ,6|. Since many 
of them make use of or want to avoid electron tunneling 
processes, their timescales are of primary importance. 

Electrons confined in semiconductor QDs interact with 
the phonon environment, which leads to loss of coherence. 
Tunneling coupling assisted by phonon degrees of free- 
dom opens an additional decoherence channel. In the case 
of doubly doped structures, the Coulomb interaction plays 
also an important role. 

In this contribution, we study phonon-assisted singlet- 
singlet relaxation in laterally coupled GaAs QDs. QDMs 
doped with two electrons are analyzed and Coulomb inter- 
action between them is included. Interaction with acoustic 



It is shown that the phonon-assisted tunneling rates reach 
much higher values (up to 160 ns^^ even at zero temper- 
ature) in comparison with other decoherence processes 
like spin-orbit coupling 0.01 ns^^). The influence of 
Coulomb interaction is discussed and its consequences 
are indicated. We calculate the relaxation rates for GaAs 
quantum dot molecules and study the dependence on 
quantum dot size, distance and offset between the con- 
stituent quantum dots. In addition the temperature de- 
pendence of the tunneling rates is analyzed. 
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phonons via deformation potential and piezoelectric cou- 
pling is taken into consideration. We show that the rates of 
phonon-assisted singlet-singlet relaxation are high (in the 
picosecond range) and strongly depend on the system pa- 
rameters and temperature. 

2 The model A QDM doped with two electrons is 
considered. The Hamiltonian of the electron subsystem reads 



2m* 



(V^ + V^)+C/(ra)+(7(rb)+y(ra,rb), (1) 



where m* = 0.07mo is the effective mass of an elec- 
tron in GaAs and F(ra,rb) = e^/(47reoei|ra — rb|) is 
the Coulomb interaction between the electrons. Here, e de- 
notes electron charge, Eq is the vacuum dielectric constant, 
and is the static relative dielectric constant. The confine- 
ment potential ?7(ra/b) for two electrons referred to as 'a' 
and 'b' is assumed to be separable: U{r) = ^m*uj'^ + 
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^m*Ljy + U{x), where the lateral potential describing 
the double QD is chosen in the form 

l/K-d/2\2 l/:c + d/2\2 

U{x) = -Uae---y^^) -{Uo + AU)e-^y^^) . (2) 

Here, d denotes the distance between QDs, a is the width of 
wave function in the x and y direction, Uo and Uq+AU are 
QDs confinement depths and the difference AU between 
them is referred to as the offset. The dynamics in the y 
and z directions is restricted to the ground states described 
by Gaussian wave functions. 

Two-particle spin-singlet states, labeled as |0) and 
are constructed from the numerically obtained single par- 
ticle wave functions, where the Coulomb interaction be- 
tween electrons is also included Q ■ 

The Hamiltonian describing the interaction of the elec- 
tron singlet states with phonons reads 



(II 



s,k 



Gsik){b 



.-fe 



H.c. 



(3) 



where k is the phonon wave vector, bl ^ and ^ are phonon 
creation and annihilation operators, respectively. Gg {k) are 
the two-electron coupling constants for different phonon 
branches s. We take into account interaction via deforma- 
tion potential and piezoelectric coupling with respective 
coupUng elements: 



DoGik), 




M,{k)g{k). 



(4) 



(5) 



2pVcsk eo£r 

Here, p denotes the crystal density, V is the normalization 
volume of the phonon modes, Cg is the speed of sound (lon- 
gitudinal or transverse), is the deformation potential 
constant for electrons and dp is the piezoelectric constant. 
The form factor Q{k) reflects the geometrical properties 
of the wave functions and Ms{k) is a function depending 
only on k orientation (for details see 13). 

We describe the properties of the phonon reservoir by 
means of the phonon spectral density 

Riuj) = (6) 

-2 |?lfe + 11^^ \Gs{k)f [S{UJ - UJs,k) + 5{l^ + LOs,k)\ , 
s.k 

where rik is the Bose distribution function. 

The phonon-assisted tunneling rates are calculated within 
the Markovian approach. To this end, we employ the Fermi's 
golden rule and obtain the rate in the form 



[AE 



(7) 



The relaxation rates are proportional to the phonon spectral 
density at the frequency corresponding to the splitting en- 
ergy AE, which is the difference between the two singlet 
states energies. 



3 Results In this section, we present and discuss the 
results for phonon-assisted tunneling rates in GaAs quan- 
tum dot molecule. In the first part, relaxation rates are cal- 
culated at r = K, at the end the temperature dependence 
is shown. 

First, the relaxation rates due to interaction with phonons 
via deformation potential are considered. The resulting tun- 
neling rates are plotted in Figs.[nAl)-(A3) for three differ- 
ent quantum dot sizes (see the sketches on top) as functions 
of the offset AU and the distance between QDs d. For 
all three configurations of QDMs, the rates are not sym- 
metric with respect to the offset (in contrast to the QDM 
doped with only one electron which results from the 
Coulomb interaction between electrons. The maximal value 
of the rate is shifted towards larger offsets with growing 
distance. For QDs close to each other, the splitting energy 
is larger than the acoustic phonon energies, which leads 
to low relaxation rates. Next, in a certain range of dis- 
tances, the energy splittings are comparable with phonon 
energies and tunneling processes are very efficient, which 
leads to the high relaxation rates. For large distances, the 
relaxation rates decrease again, which results from the de- 
creasing overlap between the two corresponding electron 
singlet wave functions. The rates are low for small offsets, 
where the Coulomb interaction prevents the electron from 
tunneling, and for large offsets, where the energy gap be- 
tween levels is very large. The values of the tunneling rate 
are smaller for larger structures since the carrier-phonon 
interaction is, in general, weaker The relaxation processes 
cover smaller parameter ranges for larger QDs since such 
structures need larger distances and the range of energies 
fitting to the acoustic phonon energies is smaller. 

In Figs.[nBl)-(B3), the tunneling rates resulting from 
piezoelectric coupling are plotted. They cover smaller pa- 
rameter ranges than rates due to deformation potential in- 
teraction since their spectral density covers only low fre- 
quencies. It is interesting to note that in the case of small 
quantum dots, this contribution is much smaller then the 
one resulting from deformation, whereas for larger QDs, 
the values of relaxation rates due to both phonon interac- 
tions are comparable. This results from the fact that the 
tunneling of an electron leads to larger change of charge 
distribution in larger structures (larger distances between 
QDs) and hence to larger piezoelectric effects. Therefore, 
this coupling is of great importance in double quantum dot 
structures in contrast to the phonon response after optical 
excitations in single QDs, where the piezoelectric coupling 
is in many cases negligible ll8ll9l [T0l . 

The total rates of tunneling involving all phonon modes 
are shown in Figs. [nCI)-(C3). For the smallest consid- 
ered QDM, the tunneling rates are high for offsets between 
AU = 17 and AU = 29 meV and for distances from 
d = 18 to d = 38 nm. They reach the maximal values of 
^ 160 ns^^ for AU ^ 23 meV and d 23 nm. For sym- 
metric QDs, the high rates cover smaller values and range 
of the offset (from AU = 14.5 to AU = 21 meV) and 
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Figure 1 Top: Schematic plots of the three considered QD sizes. (Al), (A2), and (A3) Tunneling rate assisted by phonons 
via deformation potential coupling as a function of the offset AU and distance d between QDs. (Bl), (B2), and (B3) 
Tunneling rate assisted by phonons via piezoelectric coupling. (CI), (C2), and (C3) Total tunneling rate. 



similar range but larger values of the distance (from d = 20 
to d = 40 nm). The maximum of 75 ns~^ is reached for 
AU ~ 18 meV and d ^ 26 nm. For large QDs, the high 
relaxation rates up to 25 ns^^ cover the offset range from 
AU = 7 to AU =11 meV and the distance range from 
d = 34 to d = 50 nm. 

We have shown the tunneling rates for three different 
sizes of QDs. To have a better insight in the dependence of 
the relaxation rates on the QD sizes, we present in Fig. |2] 
the rate as a function of the offset and of the lateral QD 
size a. The rates are calculated for fixed QD size in the 
growth direction, h ^ 2 nm, and the distance between QDs 
d = 30 nm at T = K. The maximal values of the tun- 
neling rate are shifted towards smaller offsets with growing 
QD size. This results from the fact that for larger structures, 
the energy differences are smaller and the Coulomb inter- 
action is weaker, so one does not need a large QD offsets 



for phonon-assisted tunneling taking place in the x direc- 
tion. 

Next, we study the relaxation rate as a function of the 
QD size in the growth direction h, which is presented in 
Fig. |3] Here, the dependence is simpler than for the dif- 
ferent QD sizes in the x direction, where the electron tun- 
neling occurs. The relaxation rates have higher values for 
smaller h and cover a wider sector of the offset. The electron- 
phonon interaction is stronger for smaller structures, since 
the frequency range of phonon spectral density is larger 
with a cutoff at ujg ~ Cs/h. 

Up to now, all the presented results were calculated at 
T = K. In Fig. |4l the tunneling rates at three different 
temperatures are shown. One can see the strong tempera- 
ture dependence: already at T = 10 K, the rates are higher 
by factor 2 than at zero temperature. This is a direct con- 
sequence of the scaling behavior of the Bose distribution 
function. 
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Figure 2 Tunneling rate as a function of the offset AU and 
the size of the QD in x and y directions a {h = 2 nm and 
d — 30 nm). 
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Figure 3 TunneHng rate as a function of the offset AU 
and the size of the QD in z direction h (a ^ 12.5 nm and 
d = 40 nm). 



4 Conclusion We have studied phonon-assisted tun- 
neling rates in semiconductor quantum dot molecules. For 
structures doped with two electrons, the singlet-singlet re- 
laxation channel was considered. The carrier-phonon inter- 
action due to deformation potential as well as piezoelectric 
coupling was studied in detail, and the importance of the 
piezoelectric effects in QDM was indicated. The depen- 
dence of tunneling rates on the QDM parameters like the 
sizes, the offset and the distance between constituent QDs, 
as well as on temperature was analyzed. 

The phonon-assisted tunneling rates reach high values, 
so comparing with the spin coherence times on a milisec- 
ond timescale 111] the process is up to several orders of 
magnitude faster The relaxation rate values are also com- 
parable with relaxation times in a single QDs jTZl . This 
shows, that one cannot neglect phonon-assisted tunneUng 
processes with such time scales. 
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Figure 4 Tunneling rate as a function of the offset AU at 
three different temperatures. 
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